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Abstract. We establish the quaternionic weighted zeta function of a graph and its Study deter¬ 
minant expressions. For a graph with quaternionic weights on arcs, we define a zeta function by 
using an infinite product which is regarded as the Euler product. This is a quaternionic extension 
of the square of the Ihara zeta function. We show that the new zeta function can be expressed as 
the exponential of a generating function and that it has two Study determinant expressions, which 
are crucial for the theory of zeta functions of graphs. 

ODD 


1 Introduction 

Zeta functions of graphs started from Ihara zeta functions of regular graphs by Ihara [8]. 
Originally, Ihara [8] presented p-adic Selberg zeta functions of discrete groups, and showed 
that its reciprocal is an explicit polynomial. Ihara also showed the logarithm of Ihara zeta 
function has an expression in the form of a generating function. Serre [13] pointed out that 
the Ihara zeta function is the zeta function of the quotient T/T (a finite regular graph) 
of the one-dimensional Bruhat-Tits building T (an infinite regular tree) associated with 
GL(2,k p ). A zeta function of a regular graph G associated with a unitary representation 
of the fundamental group of G was developed by Sunada [IB] . im Hashimoto [h] treated 
multivariable zeta functions of bipartite graphs. For a general graph, Hashimoto [BJ gave a 
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determinant expression for its Ihara zeta function by its edge matrix. Bass [2] generalized 
Ihara’s result on the Ihara zeta function of a regular graph to an irregular graph, and showed 
that its reciprocal is again a polynomial. Various proofs of Bass’ Theorem have been given by 
Stark and Terras ra, Foata and Zcilberger [5], Kotani and Sunada [9]. For the weighted type 
of zeta function of graphs, Hashimoto [2] introduced a zeta function of a graph with weights 
assigned to its edges. Stark and Terras ra defined the edge zeta function of a graph with 
weights assigned to its arcs, and gave its determinant expression by using its edge matrix. 
Mizuno and Sato HQ introduced a special version of the edge zeta function of a graph, and 
defined the weighted zeta function of a graph by using arc weights and the variable t counting 
the length of cycles. Later, they called this zeta function the first weighted zeta function in 
order to distinguish it from another zeta function which Sato defined after. It is noteworthy 
that the expressions of those zeta functions by the exponentials of generating functions were 
important to deriving the determinant expressions in the studies above. Recently, Watanabe 
and Fukumizu [[18] presented a determinant expression for the edge zeta function of a graph 
by using matrices with size of the number of its vertices. As described above, zeta functions 
of graphs have been investigated for half a century. 

On the other hand, the quaternions were discovered by Hamilton in 1843. It can be 
considered as an extension of the complex numbers. However, quaternions do not commute 
mutually in general. For many years, a number of people, for example Cayley, Study, Moore, 
Dieudonne, Dyson, Mehta, Xie, Chen, have given different definitions of determinants of 
quaternionic matrices. Detailed accounts on the determinants of quaternionic matrices can 
be found in, for example, BUS]. In order to extend the zeta function of a graph to the case 
of quaternions, we will use the approach developed by Study [T3]. The Study determinant is 
the unique, up to a real power factor, functional which satisfies the following three axioms 

P : 

(Al) d(A) = 0 <=> A is singular. 

(A2) d(AB) = d(A)d(B) for all A, Be Mat(n, HI). 

(A3) If A' is obtained form A by adding a left-multiple of a row to another row or a 
right-multiple of a column to another column, then d(A') = d(A). 

Therefore it is essential to investigate the relation between the Study determinant and the 
quaternionic zeta function of a graph. The advantage of this approach is that one can reduce 
a calculation of the Study determinant to that of the ordinary determinant. Its disadvantage 
is that the Study determinant is not an exact extension of the determinant but is rather that 
of the square of it. 

Our aim in this paper is to define the quaternionic zeta function of a graph by using 
an infinite product, namely the Euler product, and to show its essential characteristics. 
Specifically, we derive its expression in the form of the exponential of a generating function 
and determine its two types of Study determinant expressions such as in [5] (Hashimoto 
type) and in [2] (Bass type). Both types are crucial for a zeta function of a graph. In 
order to obtain these results, we explain the formal power series on a monoid algebra for 
preparations. The formal series, which includes the formal power series, is treated in detail 
in |3j Our approach in this material follows the manner in [5j ;12]. Furthermore, we give an 
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account of some relations between their exponentials and logarithms which can be found in 
Chapter IV of pE]. Up to the present, no quaternionic extension of zeta functions of graphs 
has been appeared. Hence our results can be regarded as the first successful construction of 
a quaternionic zeta function of a graph. 

The rest of the paper is organized as follows. Section 2 treats noncommutative formal 
series and monoids. The notion of Lyndon word is given, and some application of the factor¬ 
ization theorem is discussed. The conclusion (Proposition 12.3j) plays an essential role in our 
investigation. In Section 3, we define the exponential and the logarithm of noncommutative 
formal power series and give some formulas of which we make use in Section 6. In Section 
4, we explain the Study determinant of a quaternionic matrix and extend it to the matrix 
whose entries are formal power series with quaternionic coefficients. In Section 5, we provide 
a summary of the Ihara zeta function and the various zeta function of a graph, and present 
their determinant expressions. In Section 6, we define a quaternionic zeta function of a 
graph by using an infinite product, and give its expression in the form of the exponential of 
a generating function. In Section 7, we show that the quaternionic zeta functions has two 
types of Study determinant expressions. One (Theorem 17.11) is an analogy to Hashimoto 
type expression [6] and another (Theorem 17.21) is to Bass type |2j. 

2 Noncommutative formal power series and monoids 

Let R be a commutative ring with unity, and A an algebra over R. A[[fj] denotes the ring 
of formal power series in t with coefficients in A. Each element a of A[[t]] is expressed as 

a = ^Ta k t k (a k eA). 

k> 0 

A[[f]] can be equipped with the topology defined by the following manner. Let cu be the 
function defined as follows: 

uj :A[[f]]xH[[t]]—)-N U {oo} 

co(a,/3) = inf{n G N | ^ /3 n }. 

Then an ultrametric distance d u on A[[t]] is given by d LJ (a, /3) = 2~ u< ^ a '^ and a topology on 
A[[f]] is derived from d w . 

Let G be a monoid. i?[G] denotes the monoid algebra of G over R. R[G] is the set of 
formal sums z = J2 g e G z g g, where z g E R for each g E G and z g = 0 for all but finitely 
many g. The addition in i?[G] is coefficient-wise, and the elements of R commute with the 
elements of G in the multiplication. 

Let X = {xi, • • • ,Xn } be a finite nonempty totally ordered set in which elements are ar¬ 
ranged ascendingly. X* denotes the free monoid generated by X. Let < be the lexicographic 
order on A"* derived from the total order on A". For a word w = x^x^ - ■ -x lr E X*, r is called 
the length of w which is denoted by |iy|. The length of the empty word is defined to be 0. 
A nonempty word w in X* is called a Lyndon word if w is prime , namely, not a power w ,r 
of any other word w' for any r > 2, and is minimal in the cyclic rearrangements of w. We 
denote by L\ the set of Lyndon words in X*. It is well known that any nonempty word w 
can be formed uniquely as a nonincreasing sequence of Lyndon words. 
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Theorem 2.1. For any nonempty word w E X*, there exists a unique nonincreasing se¬ 
quence of Lyndon words li, l 2 , ■ ■ ■ ,l r such that w = l\U- ■ ■l r ■ 

Proof. For the proof, see for example ra- □ 

Let us consider i?[A*] [[£]]. Since (1 — /t) _1 = 1 + It + (It) 2 + • • • for every l e X* in 
R[X*] [[t]], Theorem 12.11 implies: 

> 

jf (i -/t 1 ' 1 ) -1 = wtH ’ ( 2 -f) 

ieL x wex* 


> 

in _R[X*] [[t]], where Jf means that the factors are multiplied in decreasing order. On the 
l€L x 

other hand, it follows that 


wt |m)| = {1 - (xi 4- r)t} 1 . (2.2) 

wex* 

(12.ip and (12.2j) imply the following equation: 

> 

{1 - (a?i -f-F x^t}- 1 = Jf (i - (2.3) 

ieL x 

From (12.31) . we obtain: 

Proposition 2.2. 

< 

1 - (xi H-h x N )t = ft (1 - (2.4) 

l€L x 

< 

where Jf means that the factors are multiplied in increasing order. 
i&l x 

Proof. In order to show that 

{nu-«"')-*}{ nd-# w )}=i. < 25 ) 

1&L X 1&L X 

we check that the coefficient of t r of the left hand side is equal to 1 if r = 0 and 0 if r > 0. 
Since nfeL*(! “ = Uw_l x (1 + It ^ + l 2 t 2 \ l \ H— •), the coefficient of at most r-th power 

of t in the left hand side of (12.51) is the same as that of the product: 

| ff (1 + lt w + l 2 t 2W + •••)}{ AC 1 * 

l£L x 1&L X 

|J|<r \l\<r 

for every nonnegative integer r > 0. This is a finite product since \X\ < oo and therefore is 
equal to 1. Thus (12.5|) holds. □ 
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Let [n] = {1,2, ••• ,n} with the natural order and [n] x [n] the cartesian product of 
[n] with the lexicographic order derived from the natural order on [n] . We set X = 
{x{i,j) | (i,j) E [n] x [n] } with the total order derived from [n] x [n] . For each matrix A = 
(dij) E Mat(n, A), we dehne p A to be the A-algebra homomorphism from R[X*] to Mat(n, A) 
dehned by p A (x(i,j )) = where E ,,j denotes the (i, j)-matrix unit. Let A (i,j) = a l} E ?J -. 

One can extend p A to the i?-algebra homomorphism p A from i?[X*][[f]] to Mat(n, A)[[t]] by 
defining p A (t) = t. For each Lyndon word l = (ii, Ji)(i 2 , J 2 ) e L [n]x[n] , we put 

A; = A(ii,ji)A(i 2 , j 2 )- ■ A(i r ,j r ). Then we have 

< 

In - {A(l, 1) + A(l, 2) + • • • + A(n, n)}t = (I n - A/'). (2.6) 

^-^[n] X [n] 

However we easily see that 

^ ( @XlZ2®i2i3 ‘ * — lir^irjr^iljr ^ A A+l h 0 * * *5 ^ I? 

[ O n (zero matrix of size n) otherwise, 

and that A(l, 1) + A(l, 2) + • • • + A(n, n) = A, hence we obtain: 

Proposition 2.3. Let A E Mat(n, A). Then, 

< 

I n At (In ’ '^ir—lir^irjr^XljA )> (2-7) 

(u ,jl ) • • • (*r ,jr ) e L [ n ] x [„] 
jk=ik +1 (fc=l,-,r-l) 

in Mat(n, A)[[t]]. 


3 Exponentials and Logarithms of noncommutative for¬ 
mal power series 


Let K = Q[[t]\ and X = {xi, ■ ■ ■,x ^} as in the previous section. Then K[X*] = Q[X*][[f]], 
the ring of formal power series with coefficients in Q[X*] where t is considered as a central 
element in K[X*]. Every a E Q[X*][[f]] can be written as a = Efc> o a kt k with a* E Q[X*]. 
For an a which has zero constant term, we dehne the exponential of a to be the element of 
Q[X*][[f]] as follows: 


expo = 



: a 


m\ 


(3.8) 


m> 0 


Similarly, for an a whose constant term equals 1, we dehne the logarithm of a as follows: 


logo: 


£ 

n>l 


A) 


n— 1 


(a 


n 


l) n . 


(3.9) 


If cco — 1, then a — 1 + (3 for some 0 E Q[X*][[f]] with zero constant term. Then we may 
write logcc = log(l + /l) = E n >i( — and logcc has zero constant term. Conversely, 
if 3 has zero constant term, then the constant term of exp 0 equals 1. 
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Theorem 3.1. Let a, (3 G Q[A"*][[f]]. 

(1) If «o = 1; then exp log a = a. 

(2) If (3 has zero constant term, then log exp (3 = (3. 

Proof. We shall prove (1). Letting a = 1 + (3, we have: 

1 / /_l\n 1 \ jyl 

exp log a = exp log(l + P) = Y ~~f( - Pj ■ (3-10) 

m >0 n >1 


We consider the partial sum of (13.101) as follows: 


m =0 n =1 


(- 1 )’ 


n 


-F 


(3.11) 


If k, t > r, then the coefficient of at most r-th power of (3 in (13. 10[) is equal to that in (13.111) . 
On the other hand, by the Taylor expansion we have: 


I+x — exp log(l + x) = 


(- 1 ) 


n —1 


-X 


ml 

m >0 n >1 


n 


(3.12) 


for a real number x with |cc| < 1. Similar to (13.lip , if k, t > r, then the coefficient of at most 
r-th power of x in (I3.12p is equal to that in (j3.13j) : 


k 1 i 
m =0 n =1 


(- 1 ) 


n— 1 


-X 


n 


(3.13) 


Hence the coefficient of at most r-th power of f3 in (13.lOj) is equal to that of x in (I3.12p . Since 
r is an arbitrary nonnegative integer, the coefficient of any power of /3 in (I3.10|) coincide with 
the coefficient of the same power of x in (13 .1 2|) . Moreover (3 = ^ n>1 a n t n implies that the 
coefficient of at most r-th power of t is determined by the terms of at most r-th power of f3 
in (I3.10p . Since r is arbitrary, we obtain that explog(l + (3) — 1 + /? in Q[A"*][[t]]. (2) can 
be proven in the same manner. □ 

Proposition 3.2. Assume a, (3 € Q[A"*][[t]] have zero constant terms. If af3 = (3a, then 
exp(a + (3) = (exp a) (exp (3). 

Proof. By easy calculations, we have: 

exp(a + /3) = J^ = (exp a)(exp/3). 

ra>0 k,£> 0 
k-\-i=m 


□ 


Proposition 3.3. Assume a, (3 G Q[X*][[£]] satisfy a 0 = (3 0 = 1 and put a' = a — 1, (3' = 
(3 — 1. If a' (3' = (3'a', then log (a (3) = logo + log/3. 
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Proof. Let p = log(l + a'), q = log(l + /?'). Both p and q have zero constant term. By 
Theorem 13.11 (1), expp = 1 + a', exp q = 1 + (3'. Since pq = qp, (1 + c/)(l + j3') = 
(exp p) (exp q) = exp(p + g) by Proposition 13.21 Thus log((l + c/)(l +/?')) = logexp(p + g) = 
p + q = log(l + a') + log(l + /3') by Theorem 13.11 (2). □ 

Corollary 3.4. Assume a E Q[X*][[f]] satisfies a 0 = 1 and put a' = 1 — a. Then logo -1 = 
— logo where or 1 = (1 — o ') -1 = 1 + of + (o ') 2 + • • • is in Q[X*][[£]]. 

Proof. It follows from Proposition 13.31 that: 

0 = log 1 = log((l - a')(l - o')” 1 ) 

= log((l-a')(l + a / + (a') 2 + •••)) 

= log(l — A) + log(l + OL + (oZ ) 2 + • • •) 

= log(l - a') + log((l - a') -1 ). 


□ 

Let A be a Q-algebra. We set X = {x(i,j) | ( i,j) E [n]x[n]} in the same manner as in 
the Section 2. For each matrix A = (a tJ ) E Mat(n, A), we define cr A to be the Q-algebra 
homomorphism from Q[A"*] to A defined by cr A (x(i, j)) = . It is clear that Theorem 13.11 

Proposition 13.21 Proposition 13.31 and Corollary 13.41 hold even if one replaces Q[A*] with A 
since one can extend cr A to the Q-algebra homomorphism <r A from Q[X*][[£]] to A[[£]] by 
defining cr A (t) = t. 


4 The Study determinant of a quaternionic matrix 

Let El be the set of quaternions. El is a noncommutative associative algebra over M, whose 
underlying real vector space has dimension 4 with a basis 1, i,j, k which satisfy the following 
relations: 


i 2 = j 2 = k 2 = — 1 , ij = —ji = k, jk = —kj = i, ki = —ik = j. 

For x = Xo + X\i + x^j + x%k E El, x* denotes the conjugate of x in El which is defined by 
x* = xo — X\i — X 2 j —X 3 k, and Re x = xq the real part of x. One can easily check xx* = x*x, 
(. x*) n = ( x n )*, and x~ l = x*/\x\ 2 for x 7 ^ 0. Hence, El constitutes a skew field. We call 
|x| = y/xx* = y/x*x = y/xg + x\ + x^ + x\ the norm of x. Indeed, | • | satisfies 

(1) |x| > 0, and moreover |x| = 0 x = 0, 

( 2 ) \xy\ = \x\\y\, 

(3) \x + y\ < \x\ + \y\. 

Any quaternion x can be presented by two complex numbers x = a + jb uniquely. Such 
a presentation is called the symplectic decomposition. Two complex numbers a and b are 
called the simplex part and the perplex paid of x respectively. The symplectic decomposition 
is also valid for a quaternionic matrix, namely a matrix whose entries are quaternions. 
Mat(m, n, El) denotes the set ofmxn quaternionic matrices and Mat(m, El) the set of mxm 
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quaternionic matrices. For M G Mat(m, n, H), we can write M = M s + jM p uniquely 
where M‘ S ,M P G Mat(m, n, C). M s and M p are called the simplex part and the perplex 
part of M respectively. We define if to be the map from Mat(mxn, HI) to Mat(2mx2n, C) 
as follows: 


if : Mat(mx?r, HI) —* Mat(2mx2n, C) 


M i —y 


'M s 

M p 


—M p " 

XP ’ 


where A is the complex conjugate of a complex matrix A. Then if is an M-linear map. We 
also have: 

Lemma 4.1. Let M G Mat(mxn, H) and N G Mat(nxm,H). Then 

if(MN) = ^(M)^(N). 


Proof. Let M = A + j B and N = C + j D be the symplectic decompositions of M and N. 
Then, 

MN = (A + j'B)(C + j~D) = AC + AjD + jBC + jBjB. 

Since X) = )X for every complex matrix X, we obtain: 


MN = AC - BD + j( AD + BC), 


and therefore 








AC - BD 

-AD - BC 



AD + BC 

AC - BD 


On the other hand, 







\A — B] 

rc -d] 


[AC -BD 

-AD - BC 

if(M)if(N) = 

B A 

D C 

— 

BC + AD 

-BD + AC 


Thus V>(MN) = ^(M)^(N) holds. 

Proposition 4.2. If m = n, then if is an injective M.-algebra homomorphism. 
Proof. By Lemma 14.11 if is an R-algebra homomorphism. Injectivity of if is clear. 

As in pp, one can characterize the image of if in Mat(2n, C) as follows. 
Lemma 4.3. Let J be the 2nx2n matrix defined as follows: 

7 _ r 0 -I n 

In 0 


Then, 


^(Mat(n, H)) = {Ng Mat(2n, C) | JN = NJ}, 
where N denotes the complex conjugate of N. 

Proof. For the proof, see |Ij. 


□ 


□ 


□ 




















In [15] , Study defined a determinant of a nxn quaternionic matrix which we denote by 
Sdet as follows: 

Sdet(M) = det(^(M)), 

where det is the ordinary determinant. Sdet is called the Study determinant. 

Sdet can be extended to the one for matrices whose entries are formal power series with 
coefficients in El in the following way. Let H[[£]] be the ring of formal power series with 
coefficients in EL We notice that t is a commuting indeterminate for H, that is, th = ht for 
any h G EL For a = Yhk>o a ^ k EI[[t]], the conjugate of a is defined by a* = J2k>o a l^ k - 
Furthermore, let ag. = af. +joi p be the symplectic decomposition of a*.. Then the symplectic 
decomposition of a is defined by 


a = a s + ja p = a%t k + j a p t k . 

k> o k> o 


One can extend -0 to the injective homomorphism ip t of M-algebras from Mat(n, H)[[f]] to 
Mat(2n, C)[[£]] by defining = t. On the other hand, since Mat(2n, <C)[[i]] = Mat(2n, C[[t]]), 
one can also extend the determinant det : Mat(2n, C)—>C to det* : Mat(2n, C)[[f]]—»C[[i]] 
in the natural manner since det is a polynomial of entries of a matrix. Both fy and det* are 
continuous, thus 

det: Mat(n, H)[[t]]—»C[[i]] 

is also continuous with respect to the topology of formal power series. We call det* the 
Study determinant for Mat(n, H[[f]]) and denote by Sdet*. We notice that the restriction of 
Sdet* to Mat(n, El) yields Sdet. 

Before stating properties of Sdet*, we mention some useful facts. 

Lemma 4.4. If A, B,C,D are complex square matrices with same size. Suppose that A is 
invertible and AC = CA, then 


det 


A B 
C D 


det (AD - CB). 


Proof. For the proof, see for example 


□ 


Remark 4.5. We notice that the Lemma f.f can be extended to matrices whose entries 
belong to C[[t]] immediately. 

Lemma 4.6. 

, 0t(Mat(n, H[[*]])) = {N G Mat(2n, C[[t}}) | JN = NJ}. 

Proof. By Lemma 14.31 it follows that 

3M = M kt k G Mat(n, H [[t]]), 0 t (M) = N = N k t k 

k>0 k>0 

-v4 ^(Mfc) = Nfc for \/k >0 JNfc = N^.J for \/k >0 
JN = NJ. 


Hence, the assertion holds. 


□ 


9 





Lemma 4.7. Let a = J2 k>0 &kt k £ H[[t]], and a k = of + ja p the symplectic decomposition. 
Then, 

aa* = a*a = ( a k a e + a k a e)t m = ®- S ot s + a p a p G C[[£]]. 

m> 0 k+l=m 

Proof. By a direct calculation, we have: 

aa * = (^2( a k + j a k) tk )(52( a e + jafYt*) 

k> 0 ^>0 

= (J2(^k + j a k) tk )C^2( a i - a ij)t e ) 

k> 0 ^>0 

m>0 k+t=m 

= ^2 ( a k a £ + a k a f) tm 

m> 0 k+i=m 

= a s a s + a p a p . 

Furthermore, a s a s + a p a p belongs to C[[£]]. Similarly, the same conclusion holds for a* a. 

□ 

Lemma 14.71 immediately implies: 

Corollary 4.8. For all a,/3 £ H[[£]], aa*/3/3* = (3/3*aa*. 

Proposition 4.9. (i) Sdetf(M) G 1 R[[t]\ for M G Mat(n, EI[[t]]) 

(ii) The constant term of Sdett(M) is equal to 0 <=>• M is singular. 

(hi) Sdett(MN) = Sdet t (M) Sdett(N) for M, N G Mat(n, B[[t]]). 

(iv) If N is obtained from M by adding a left-multiple of a row to another row or a right- 
multiple of a column to another column, then Sdett(N) = Sdet*(M). 

(v) If N is obtained from M by exchanging a row with another row or a column with 
another column, then Sdet t (N) = Sdet t (M). 

(vi) Sdett(aM) = SdeL(Ma) = (aa*) n Sdett(M) for M G Mat(n, lH[[f]]) and a G H[[f]]. 

(vii) //Mg Mat(n, H[[t]]) is of the form: 



Ai 

* 

* 


'Ai 

0 • 

• 0" 

M = 

0 

A2 

* 

or 

* 

A2 

0 


0 

0 • 

V . 

3 


* 

* 

1 

g 


Then, Sdett(M) = Uti W- 
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(viii) Let A be an m x n matrix and B an n x m matrix. Then 

Sdet t (I m - AB) = Sdet t (I n - BA). 
Proof, (i) It follows from Lemma [4.61 that 


dett(N) = det t (J -1 NJ) = det t (N) = det*(N), 

for any N G HI[[t]])). Thus Sdet t (M) must belongs to M[[t]]. 

(ii) It is known that a matrix over a commutative ring is invertible if and only if its 
determinant is invertible, and that a formal power series is invertible if and only if its constant 
term is invertible. Hence, for any M G Mat(n,H[[t]]), the constant term of detf^t(M) is 
equal to 0 44- detf , 0t(M) is not invertible ^ t (M) is singular. Taking conjugate and 
inverse of i/y (M) J = J'0 t (M), we have: 

J^(M)" 1 = @pJ. 

Hence by Lamina 14.61 ^(M) -1 G ^*(Mat(n, HI[[t]])) and it follows that ^(M) is singular 
<G> M is singular. 

(iii) Since A is an algebra homomorphism, (iii) immediately holds. 

(iv) Let b = bi + jb 2 G H[[i]] (bi,b 2 G C[[t]]). It suffice to show that for k ^ i, the nxn 
matrix B ke (b) = l n + bFj ki satisfies Sdet t (B^(6)) = 1. Since four blocks of 


A(B M {b)) = 


I n + — b 2 ^M 

b 2 ^k£ In + b\Ek£ 


commute with each other, we can apply Remark 14.51 and obtain: 

detf^(Bfc£(6)) = dett((I n + 6iE^)(I„ + 6iE^) + b^^b-zEke) — 1- 

(v) We prove for row exchange. Putting N = PM for the permutation matrix P corre¬ 
sponding to some transposition, it follows that: 


Sdett(N) = Sdetf(PM) = dep 


P O 
O P 


Am 


= deti 
= Sdeti(M). 


P O 
O P 


dett('0 t (M)) = det(P) 2 Sdetj(M) 


(vi) Let M G Mat(n, H[[£]]) and a = a s + ja p G HI [[t]\ ( a s ,a p G C [[t]}) be the symplectic 
decomposition. Then by using Lemma [4.11 Remark 14.51 and Lemma [4.71 we have: 

Sdet t (aM) = det t (A(oM)) = det t (ip t (cd n )A (M)) 


= det* 


« 5 In 

a p L 


-Of I r 
a s l n 


det t (Am) 


= dett((a' s 'l n )(a 5 In) + (a p I„)(a p I„))Sdet t (M) 
= det^c^a: 5 + a p a p )I„)Sdet t (M) 

= (acP) n Sdet t (M). 
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In the same way, we can deduce Sdet^Ma) = (aa*) n Sdet t (M). 

(vii) For a 2nx2n matrix N and any two subsets / = {ii, i 2 , • • •, i r }, J = {ju j‘ 2 , • • -, .7. s } 
of [2 n], N /J denotes the submatrix obtained from N by deleting ii, i 2 , ■ ■ •, i r th rows and 
• • •, j s th columns. Then by definitions of Sdet t and fy, we get the following: 


Sdett(M) 


deti(^(M)) 


deti 


'M s 

M p 


—M p 


Af * 
0 Af 


* —Af * 

* 0 — Af 


* 

* 


deti 


0 0 
Af * 
0 Af 


Af 0 0 

* Af _*_ 

* 0 Af 



* 


* 


[0 0 ••• Af 0 0 ••• Af J 

Af^ t (M) {1}{1} + (-l) n+2 Af^(M) {n+1}{1} 

Af Af^(M) {1 ’ n+1}{1 ’” +1} + Af Af^(M) {1 ’ n+1}{1 ’ n+1} 
(Af Af + Af Af)^(M) {1 ’ n+1}{1 ’ n+1} 
AiA^t(M) {1 ’ n+1}{1 ’ n+1} 
AiA*A 2 A;^(M) {1 ’ 2 ’ n+1 ’ n+2}{1 ’ 2 ’ n+1 ’ n+2} = • • • 


nw 

i =1 


(viii) Let A be an m x n matrix and B an n x m matrix. Then, by the definition of the 
Study determinant, we have: 

Sdet t (I m - AB) = det f (r/y(I m - AB)) = det t (I 2m - r/y(A)?/y(B)) 

= det t (I 2 „ - V't(B)'0t(A)) = deti('0t(I n - BA)) 

= Sdeti(I n — BA). 


□ 

Remark 4.10. Sdet t is not multilinear as det t is. Furthermore, Sdet t ( p M) = Sdet t (M) 
does not hold in general where is the transpose of M. 

5 The Ihara zeta function of a graph 

Let G = (I/(G), E(G )) be a finite connected graph with the set V(G ) of vertices and the set 
E(G) of unoriented edges uv joining two vertices u and v. We assume that G has neither 
loops nor multiple edges throughout. For uv G E(G), an arc ( u , v) is the oriented edge from u 
to v. Let D{G) = { (u,v), (v,u) \ uv G E(G)} and |R(G)| = n, |_E(Cr)| = m, \D(G)\ = 2m. 


12 



For e = (u, v) G D(G), o(e) = u denotes the origin and t(e) = v the terminal of e respectively. 
Furthermore, let e _1 = (v,u) be the inverse of e = (u,v). The degree degu = degc v of a 
vertex v of G is the number of edges incident to v. For a natural number k, a graph G is 
called k-regular if deg <3 v = k for each vertex v of G. A path P of length £ in G is a sequence 
P = (ei, • • • , ef) of £ arcs such that e* G D{G) and tief) = o(e;+i) for i G {1, • • • , £ — 1}. 
We set o(P ) = o(e i) and t(P) = t(ef). |P| denotes the length of P. We say that a path 
P = (ei, ■ • ■ ,ef) has a backtracking if e i+1 = ef 1 for some i(l<i<£ — 1). A path P is said to 
be a cycle if t(P) = o(P). The inverse of a path P = (ei, ■ • • , ef) is the path (ef 1 , • • • , ef 1 ) 
and is denoted by P _1 . 

An equivalence relation between cycles is given as follows. Two cycles C\ — (ei, • • • , ef) 
and C 2 — {fi, , fi) are said to be equivalent if there exists k such that fj = 6 j + k for all 
j where indices are treated modulo £. Let \C] be the equivalence class which contains the 
cycle C. Let B r be the cycle obtained by going r times around a cycle B. Such a cycle is 
called a power of B. A cycle C is said to be reduced if both C and C 2 have no backtracking. 
Furthermore, a cycle C is said to be prime if it is not a power of a strictly smaller cycle. Note 
that each equivalence class of prime, reduced cycles of a graph G corresponds to a unique 
conjugacy class of the fundamental group 7 t 1 (G, v) of G at a vertex v of G. 

The Ihara zeta function of a graph G is a function of t G C with \t\ sufficiently small, 
defined by 


Z(G,t) = Z G (t) 


(1 


where [C] runs over all equivalence classes of prime, reduced cycles of G. 

Let B = (B e /)e,/eD(G) and Jo = (J e f)e,feD(G) be 2 m x 2 m matrices defined as follows: 


Be/ - 


1 if t(e) = o(f), J = fl if / = e \ 
0 otherwise, 0 otherwise. 


Then the matrix B — Jo is called the edge matrix of G. 

Theorem 5.1 (Hashimoto[B]; Bass[2]). LetG be a finite connected graph. Then the reciprocal 
of the Ihara zeta function of G is given by 


Z(G, t)- 1 = det(I - t(B - J 0 )) = (1 - t 2 ) r ~ l det(I — tA + t 2 (D - I)), 


where r and A are the Betti number and the adjacency matrix of G, respectively, and D = 
(dij) is the diagonal matrix with da = degr; i; V(G) = {vi, ■ ■ ■ ,v n }. Z (G,t) has an expression 
as the exponential of a generating function as follows: 

Z(G,«) = exp(^^V), 

k> 1 

where Nk is the number of equivalence classes of reduced cycles of length k. 

The matrix D is called the degree matrix of G. 

A weighted zeta function of a graph was first defined in Hashimoto [Bj by giving a weight 
for each edge of a graph. After that, it was generalized in Stark and Terras ra by giving 
a weight for each arc of a graph as follows. Let D(G) = {e\, • • • , e m , e m+ i, • • • , e 2 m } where 
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e m+ i = ef 1 for 1 < i < m. For each arc e* (« = 1, • • • , 2m), associate a complex number rq, 
and set u = (tq, • • • , u 2m ). Let g(C) — u n ■ ■ ■ u ik for a cycle C = (e il} • • • , e ik ). Then the 
edge zeta function Cg( u ) °f G is dehned by 


Cg(u) 


(i -g{C)r\ 


where [ C} runs over all equivalent classes of prime, reduced cycles of G. 

Theorem 5.2 (Stark and Terras[T4]). Let G be a connected graph with m edges. Then 
Cg(u )- 1 = det(I 2m - (B - J 0 )U) = det (I 2m - U(B - J 0 )), 


where U = diag(rq, • • • ,U 2 m ) is the diagonal matrix. 

In (TTj, the multivariable analogue N k ( u) of N k was dehned and used to obtain the 
determinant expression. Mizuno and Sato HU dehned the weighted zeta function of a graph 
G, and gave a determinant expression by the weighted matrix of G. Later, they called this 
zeta function the hrst weighted zeta function in order to distinguish it from another zeta 
function which Sato dehned after. Assigning a complex number w(e) to each arc e G D(G), 
we dehne W = ( w uv ) u ,vev(G) to be the n x n matrix as follows: 


f w(e) if e = (u,v) G D(G), 
\ 0 otherwise, 


W is called the weighted matrix of G. For each path P = (e^, e.; 2 , • • •, e^), the norm w(P) 
of P is dehned by w(P) = w(ei 1 )w(ei 1 )- ■ ■w(ei d ). The (first) weighted zeta function of G is 
dehned by 


Z (G,w,t) 


(1 -w(C)tM)-\ 


c 


(5.14) 


where [C] runs over all equivalent classes of prime, reduced cycles of G. 


Theorem 5.3 (Mizuno and SatofTT]). Let G be a connected graph with n vertices. Suppose 
that w(e~ l ) = w(e)~ 1 for e G D(G). Then 

Z (G, w, t)- 1 = (1 - t 2 ) m ~ n det(I„ - tW + t 2 (D - I n )). 


Z (G,w,t) has an expression as the exponential of a generating function as follows: 

Z (G, w, t) = exp ( ^ - w ( c ') s t |C|s ) i 

[G] s>l S 

where [C] runs over all equivalent classes of prime, reduced cycles of G. 

Later, Watanabe and Fukumizu [TS] gave a determinant expression for the edge zeta 
function of a graph with no additional condition w(e~ 1 ) = w(e)~ l . Rearranging arcs so that 
D(G) = {ei, ef 1 , • • • , e m , e~f}, and putting Ui = w(ef), U turns out to be 

w(e i) 0 

w^f 1 ) 
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In [18], two nxn matrices A = (A uv )u, v eV(G) and D = {JH av ) u ,v&v(G) were defined as follows: 


A 

-f'-uv 


D 


UV 


1 

0 


w(u, v ) 

w(u, v)w(v, u ) 


if (u, v ) G D(G), 
otherwise, 


v- w(e)w{e *) 
uv ^ 1 — w(e)w(e~ 1 )' 

eeD{G) v y V ) 

o(e)=u 


where 6 UV denotes the Kronecker’s delta. They derived their determinant expression with A 
and D instead of the weighted matrix and the degree matrix of G. 

Theorem 5.4 (Watanabe and Fnk nmi zn|18j). 

m 

Cg(u)- 1 = det(I 2m - U(B - Jo)) = det(I n + D - A) JJ(1 - w^w^ef 1 )) 

i =1 


6 The quaternionic weighted zeta function of a graph 

We continue with the notations in former sections. We shall define the quaternionic weighted 
zeta function of a graph. First we extend w(e) to a quaternion. Let w be a map from 
V(G)xV(G) to HI such that w(u,v ) = 0 if (u,v) ^ D(G). We may write w(e) instead of 
w(u,v ) if e = (u,v) G D(G). We call w(e) the quaternionic weight of e G D{G). For a 
path P = (ei, • • ■ , eg) of G, the norm (or the quaternionic weight ) w(P) of P is defined by 
w(P) = w(ei)w(e 2 )- ■ ■w(eg). We arrange the arcs e l7 e 2 , ■ ■ -, e 2m such that e m+ k = e ^ 1 for 
k — 1, ■ • •, m. Now, we define the quaternionic weighted zeta function of G to be the element 
of H[[t]] as follows: 

Z m (G,w,t) = Jj|(l-w(C , )t |c ' l )(l-u;(C')t |c ' l )*| , (6.15) 

c 

where C runs over all reduced cycles C = e^e^- • -e^ such that i\i 2 - ■ 4 r GL[ 2m ]. If w(e) — 1 
for any e G D(G), then the quaternionic weighted zeta function of G is reduced to the square 
of the Ihara zeta function of G. 

Remark 6.1. (i) The primeness is not required explicitly since i\i 2 - ■ -i r GL[ 2m ] implies the 

primeness of e tl e 82 ■ • -e lr . 

(ii) (16.15[) does not depend on the order in which inverses of (1 — ta(C)tl c 'l)(l — w(C)t^)* 
are multiplied since (1 — ta(C')fl c 'l)(l — u’(C)fl c 'l)* G C[[t]] by LemmaIndeed, the 
inverses are in M[[t]] since (1 — ta(C)fl c 'l)(l — w(C)t^ c ^)* = 1 — 2 {Reta(C')}tl c l + 
\w{C)\H 2 \ c \. 

We shall show that Z w (G,w,t) can be expressed as the exponential of a generating 
function. We recall that Theorem 13.11 Proposition 13.21 Proposition 13.31 and Corollary 13.41 
hold even if we replace Q[A*][[t]] with H[[£]]. By Remark 16.11 we have: 

Z m (G,w,t) = ]^[ |l - 2{Reu'(C)}f lc| + |w(C)| 2 t 2|G| | . 

c 
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Since ac = 2{Re w{C)}t^ — \ w(C)\ 2 t 2 \ c \ is an element of M[[£]] and has zero constant term, 


{1 - 2{Rew(C')}t |c| + HC)| 2 t 2|c| }- 1 = {1 - ac}' 1 = 1 + X«c 

k> 1 

is also an element of M[[f]] and fie = Ylk >i a c has zero constant term. Since the constant 
term of Z H (G, w,t) is equal to 1 , we can derive the following: 


logZ m (G,w,t) = log na+fh) = X log ( 1+ ^) 

c c 

= X 1 ° s (( 1 - «c) _1 ) = - X 1 ° s ( 1 “ a c) 

c c 

= -X log {( 1 _ w(C)t lcl )(l - w(C')*t |c| )} 
c 

= — X { l°g(l — w(C)t^) + log(l — w(C)*t ^)| 

c 

= ~EE —{“’(C)" + (“’(C)-)"}«"' C| 

C n> 1 


EE 


2Re(w(C)") J „ lcl _ 


n 


< n|c| = E E 


2Re(«.(C)") J „| C | 


C n>1 n >1 C 

Therefore log Ze(G, w, t) has zero constant term, and it follows that: 


n 


Theorem 6.2. Z a (G,w,t) has the expression as the exponential of a generating function as 
follows: 

2Re(tt,(CT) ; „| C r 
n 


Z h (G, w, t ) = exp (EE 


n> 1 C 


where C runs over all reduced cycles C = e^e^- • -e ir suc/i that iii 2 - ■ -vG^m]- 


We can also express Z^(G,w,t) as the square of the exponential of a generating function 
as follows: 

Re (w(C) w ) r | C |N ^ 2 
n 


Zh (G, w, t) — | exp (EE 


n> 1 C 


7 Study determinant expressions for Z u(G,w,t) 

Consider an n x n quaternionic matrix W = (W U v)u,v&v(G) G Mat(n, H) with (u, u)-entry 
equals 0 if (u,v) ^ D(G). We call W a quaternionic weighted matrix of G. Furthermore, let 
w(u,v) = W uv for u,v G V(G) and w(e) = w(u,v ) if e = (u,v) G D(G). 

For a quaternionic weighted matrix W of G, we define two 2 m x 2m matrices B u , = 
(By) eJeD(G) and = (jJ } ) e ,/ e r>(G) as follows: 

B W = f ^(e) if t(e) = o(f), y™) = f w(e) if / = e -1 , 

0 otherwise, [0 otherwise. 
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We call B^ — 3 W the quaternionic weighted edge matrix. Then the Study determinant ex¬ 
pression of Z h (G, w, t) of Hashimoto type by use of the quaternionic weighted edge matrix 
is stated as follows: 


Theorem 7.1. 

= Sdet t (I 2m - (B^ - J w )t). 

Beforehand, we mention that the right hand side does not depend on the way of arrang¬ 
ing arcs which determines row and column indexes. Indeed, if we rearrange arcs, then it 
turns into Sdet t (P(I 2m — (B m — J„,)t)P _1 ) for some permutation matrix P. However, from 
Proposition 14.91 (v). this is equal to Sdet t (I 2m — (B m — 3 w )t). 

Proof. Set A — H in Proposition 12.31 Then we can take Study determinants of both sides 
in (12.71) . so that 


Sdett(I 2m - At) 


Sdetj ^ | | (In 


(u ,31) ■ ■ ■ (V Or)&L [2m] x [2m] 
jk=ik+i (fc=l,—,r-l) 


(7.16) 


n 


Sdet t (I n ®U *2 ®* 2*3 ’ ’ )' 


(il ,jl ) — (V ,jr) £L [ 2m ] x [2m] 

jk=ik+1 (fc=l,—,1 1) 

We notice that the last formula does not depend on the order in which Sdet t are multiplied 
since Sdetf take values in M[[t]]. It follows from Proposition 14.91 (vii) that if j r = i±, then 


Sdet t (I n ®n* 2 ®i 2 * 3 ’ ' ) 

= (1 — ai x i 2 ai 2 i 3 - ■ -ciiriA )(1 — - a ili2 ai 2 i 3 - ■ -CLiriA ) j 


and otherwise, 


Sdetain (^iii 2 (^i 2 i3' ' '^trjr Ip 1 jr ( ) T 


Putting A = B^ — 3 W , then 


Q*ij ^eiej 


w(e.j) if t(ei) = o(ej ) and e^e i 1 , 
0 otherwise. 


Therefore, (17.161) yields: 

Sdeti(I 2m (B m J«;)t) 

= Sdett(ln~w(ei 1 )w(ei 2 )---w(e ir .)E ilil t r ) 

(U,i 2 )-”(*r,u)ei[ 2 m]x[ 2 m] 

£i 2 '" e i r ireduced cycle 

= (l-w(e il )---w(e ir )t r )(l-w(e il )---w(ei r )t r )*. 

(U,i2)”-(*r,n)ei[2m]x[2m] 
ei 1 ei 2 •••€i r ireduced cycle 
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Each Lyndon word (ii,i 2 )- • in L[ 2m ] x [ 2 m] corresponds to a Lyndon word Ri 2 - ■ ■i r in 

L[ 2 m] bijectively. Hence we have: 

Sdeti(I 2m (B u J w'jt') 

= JJ (1 — w^e^)- ■ ■w(e ir )t r )(l — w^e^)- ■ ■w(e ir )t r )*. 

ili2-'ir&L[2m] 
ei 2 ”'^i r :reduced cycle 

Remark 16.11 (ii) implies the inverse of the right hand side is Z ^(G,w,t). Hence we obtain 
the conclusion as desired. □ 


Finally, we shall show the Study determinant expression of Bass type. We define n x n 
matrices W = (W U v)u,vev(G) and D = {B U v)u,vGV{G) to be as follows: 


W 


uv 


D 


UV 


J (1 — w(e)w(e 1 )t 2 ) 1 w(e ) if e = (u,v)eD{G) 
1 0 otherwise, 

S uv ^2 (1 - w(e)w(e~ 1 )t 2 )~ 1 w(e)w(e~ 1 ). 


e&D(G) 

o(e)=u 

Then the Study determinant expression of Bass type is stated as follows: 

Theorem 7.2. 


Z a(G,w,t) 1 

m 

= Sdet t (I„ - tW + t 2 D) ]^[(1 - w(e i )w(e^ 1 )t 2 )( 1 - w(e i )w(eY 1 )t 2 )*. 

1=1 


Proof. We readily see that: 


Sdet^(l2 m t(T3 w J^)) Sdett(l2 m t3 w t’^uf) 

= Sdet t (I 2m - tB w (I 2 m + iJu,) -1 ) 
x Sdet t (I 2m + tJ w ). 


Now, let K 
as follows: 


(K ev ) eeD(G)-,vev(G) and L = {B ev ) e&D(G) . v&v(G) be two 2 m x n matrices defined 


J w(e) if t(e) = v, 
) 0 otherwise, 


f 1 if o(e) = v, 
[ 0 otherwise. 


Then one can easily check that: 


K t L = B„,. 


Thus, from Proposition 14.91 (viii). we get the following: 


Sdet t (I 2m - t(B w - J w )) = Sdet t (I 2m - tK T L(I 2m + tJ w ) *) 

x Sdett(I 2m + tJ w ) 

= Sdet((I„ - t T L(I 2m + tJ«,) _1 K) 
x Sdet t (I 2m T tJ w ). 


(7.17) 
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Arranging arcs so that D(G) = {ei,e 1 1 , ■ ■ ■ , e m ,e m 1 }, we have from Proposition 14.91 (iv), 
(vii) that: 


Sdet f (I 2m + LJJ = Sdett 


= SdeL 


1 w(e\)t 0 0 

w(ei l )t 1 0 0 

0 0 1 w(e 2 )t 0 

0 0 w(e^ 1 )t 1 0 

; : : o 

1 — w(ei)w(e^ 1 )t 2 w(ei)t 0 0 

0 1 0 0 

0 0 1 — w(e 2 )w(e 2 1 )t 2 w(e 2 )t 0 

0 0 0 1 0 

: : : 0 

-1\ .2\* 


IK 1 - w(ei)w(ei ^(l - w(ei)w(ei x ) t 2 ) 


1=1 


Let X(e) be the 2x2 matrix defined as follows: 

1 w(e)t 


X(e) — 

Then, the inverse of X(e) is given by 


w(e 1 )t 1 


(7.18) 


X(e)- 1 


Therefore, 


(1 — w(e)w(e 1 )t 2 ) 1 
— (1 — w(e~ 1 )w(e)t 2 )~ 1 w(e~ 1 )t 


— (1 — w(e)w(e *) t 2 ) 1 w(e)t 
(1 — ^( e -1 )^^)^ 2 ) -1 


(I 


2m 


+ tJ w ) 


'X(ei)- 1 O ■■■ O' 
O X(e 2 )" 1 

O 

. o ••• O X(e m ) -1 _ 


If e = (u,v) E D(G), then we have by a calculation that: 


( T L(l2 m + tJw) 1 K) U „ — ( T L) u j((l 2 m + tj w ) 1 )jj/K f'v 

f,feD(G) 

— ( T L)ue((l2m + tJ w ) 1 ) eg K CT 
= (1 — w(e)w(e~ 1 )t 2 )~ 1 w(e). 

If u = v, then we have: 

( T L(l2m + tj w ) 1 K ) uu — yy ( T L)u/((l2m + tJ w ) ^ //'K/'u 

f,f'eD(G) 

— ( T L) lie ((l2m + tj w ) 1 ) ee _ 1 K e -l M 

e&D(G) 

o(e)=u 

= — yy (1 — w(e)w(e~ 1 )t 2 )~ 1 w(e)w(e~ l )t. 

ee D(G) 
o(e)=u 
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Hence, it follows that: 


(7.19) 


T L(I 2m + tJJ-'K = W - tD. 

Combining (17.17)) . (17.18)) and (17.19)) , we obtain: 

Sdet t (I 2m t(B w Jtt;)) 

m 

= Sdet t (I„ - tW + t 2 t>) JJ(1 - w{e i )w(ej 1 )t 2 ){l - w{e i )w(ej 1 )t 2 )*. 

2=1 

Now, the assertion follows from Theorem 17.11 □ 
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